Quantum Creep in Yi x Pr x Ba 2 Cu 3 07 <$ Crystals: Magnetic Relaxation and Transport. 
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We report transport and magnetic relaxation measurements in the mixed state of strongly underdoped 
Yi- x Pr x Ba2Cus07-s crystals. A transition from thermally activated flux creep to temperature independent 
quantum flux creep is observed in both transport and magnetic relaxation at temperatures T < 5 K . Flux trans- 
former measurements indicate that the crossover to quantum creep is preceded by a coupling transition. Based on 
these observations we argue that below the coupling transition the current is confined within a very narrow layer 
beneath the current contacts. 



Quantum flux creep in superconductors is the only ex- 
perimentally accessible phenomenon in which a macro- 
scopic metastable system, such as a persistent supercur- 
rent, relaxes coherently without thermal activation. In 
all other metastable systems, the relaxation proceeds as 
a sequence of a large number of uncorrelated microscopic 
steps, requiring thermal activation over an energy bar- 
rier. Yet, the evidence of magnetic relaxation and re- 
sistance that do not extrapolate to zero at a tempera- 
ture T — > has not reached a point where experimental 
data begin to form a cohesive picture of the phenomenon. 
While the non-vanishing magnetic relaxation has been 
observed both in single crystals and thin films IT] |3] , the 
non-vanishing temperature independent resistance had 
been observed only in ultrathin films and multilayers 
[0-p|. This has contributed to the assertion that the 
temperature independent resistance in films and the non- 
vanishing low-temperature magnetic relaxation in single 
crystals are unrelated phenomena. 

A second question, closely related to quantum creep 
in high-T c superconductors, is the coupling transition in 
layered crystals. Several groups have arrived at conflict- 
ing conclusions. Safar et al. M reported that the voltages 
generated by the motion of vortices on the opposite faces 
of the YBa-iCu-zOi-s sample (flux transformer method) 
converge. However, such a strong manifestation of the 
coupling transition appears to be the exception rather 
than the rule. Other groups have observed that the volt- 
ages generated on the opposite faces of Bi 2 Sr2CaCu20s 
[p|-|lli an d YBaiCu?,07-£, O] crystals diverge with low- 
ering temperature, rather than converge. 

Here we report the results of both magnetic relaxation 
and transport measurements on strongly underdoped 
crystals of Y 0A7 Pr . 5S Ba 2 Cu z O 7 -s (T c « 17 - 21 A'). 
Our transport data indicate that the onset of quantum 
creep is preceded by a a coupling transition which leads 
to non-ohmic dissipation. The transition to temperature 
independent creep takes place at lower T and only in a 
very thin layer below the current contacts; the rest of the 



sample continues to exhibit thermally activated creep. A 
picture of the low temperature dissipation of the trans- 
port current that arises from these observations is that 
the sample is divided into two macroscopic regions: a 
very thin layer, a few unit cells thick, near the upper 
face, where the current contacts are located, which car- 
ries most of the current, and the rest of the sample which 
remains mostly undisturbed by the current. These two 
macroscopic regions are decoupled from each other. In- 
side these layers the vortices are coherent, with the corre- 
lation length comparable to the thickness of the respec- 
tive layer. Magnetic relaxation measurements confirm 
the transition to T-independent relaxation at approxi- 
mately the same temperatures as in transport. 

Two strongly underdoped twinned single crystals of 
YoA7Pfo.53Ba2Cu307-s with T c rs 17 and 21 K, re- 
spectively, were prepared as described in Ref. |Tj|. The 
first sample was used for transport measurements with 
the "flux transformer" contact configuration [inset to 
Fig. H(a)]: The current I was injected through the con- 
tacts on one face of the sample and the voltage drop 
between the voltage contacts on the same (primary volt- 
age V p ) and the opposite (secondary voltage V s ) faces 
was measured for temperature, total current, and mag- 
netic field H , applied parallel to the c-axis, in the ranges 
1.9 if < T < 25 K, 0.3 fxA < I < 2 mA, and 
0.2 T < H < 9 T. Magnetic relaxation measurements 
were performed on the second crystal using a SQUID 
magnetometer. The crystal was cooled in zero field; a 
field H + AH (AH = 0.3 T for all H) was applied par- 
allel to its c-axis, and then the field was reduced to H . 
The decay of the resultant paramagnetic moment was 
monitored for several hours (« 10 4 s) in constant H . 

Figure [j](a) gives an overview of the T and H depen- 
dence of V p and V s normalized to I and H . The con- 
vergence of these curves at T rs 9 K indicates a regime 
of free flow of vortices, where V p , s /I oc H. In a region 
below 9 K, V p and V s exhibit activated T— dependences 
with field-dependent activation energies. At even lower 



10 4 



10 J r 



g 


10 2 


a 






m 1 


X 

* 





10° r 



10' r 



■ i i i I i i i i I i i i i I i i i i I i i i i I i i i i 

unpinned 
flux flow 



12 3 4 

11 L - 




-, 10 4 



10 



w 



-* 10' 



-, 10 u 



^ Quantum 

• Creep 




10 3 



3. < 






10 



0.2 0.3 0.4 0.5 0.6 
1/T(K" 1 ) 

FIG. 1. (a) Primary V p and secondary V s voltages, nor- 
malized to the total current / and magnetic field H , plotted 
vs 1/T for five magnetic fields (0.2, 0.4, 0.6, 0.8 and 1 T). The 
slope decreases with increasing field. Inset: Contact con- 
figuration used in the measurements, (b) Decay time tj of 
magnetic moment (see definition in text) vs 1/T for several 
values of magnetic field (H = 0.1, 0.2, 0.6, 0.8, and 1.2 T). 
The slope decreases with increasing field. The straight line 
extrapolations of the Arrhenius-type dependence converge at 
T cr « 18 — 19 K and r esc ~ 1 — 2 s. The saturation of Td at 
the level 10 11 — 10 12 s is due to quantum creep. 

T, the primary signal becomes T-independent and scales 
approximately with H , while the secondary voltage re- 
mains thermally activated. 

Figures |2|(a) and |2|(b) represent an expanded view of 
V p>a vs 1/T for H = 0.6 T, I = 100 fiA and H = 1.5 T, 
J = 10 fiA, respectively. The activation energies E p>s = 
—din V PlS /d(l/T) of both V p and V s change suddenly at a 
temperature T* (H). In addition, the dissipation becomes 
strongly non-ohmic below T*. As the inset to Fig. g(a) 
shows, the primary resistance R P (T) = V p /I is current 
independent only above T*. The activation energies E p<s , 
which are equal and current-independent at T > T* , be- 
come current-dependent at T < T* . At low current, the 
activation energy below T* is greater than above T* 
[Fig. 0(b)] so that the curves V PiS (l/T) have a downward 
curvature. However, the activation energy decreases with 
increasing current and E p (I) becomes smaller than it is 
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FIG. 2. Arrhenius plots of the primary V p and secondary 
V s voltages measured in two different fields and currents: (a) 
H = 0.6 T and I = 100 fiA, and (b) H = 1.5 T and 
I — 10 fiA. Inset to (a): Expanded view of the primary 
resistance V p /I for two values of the total current (lfiA and 
250fiA) in an applied magnetic field H — 0.2 T. Inset to (b): 
Schematic representation of the current distribution. The cur- 
rent flows mainly in the upper (shaded) layer. The vortices 
are coupled within each layer, but the layers are decoupled 
from each other. 

at T > T* , so that V p (l/T) acquires an upward curvature 
[V p in Fig. ||a) an d the inset]. 

This evidence indicates that the vortices undergo a 
coupling transition at T*. At T > T* , the dissipation 
mechanism is activated hopping of 2D pancakes over po- 
tential barriers since the activation energies are the same 
for the primary and secondary voltages, in spite of the 
nonuniform current distribution. When the pancakes be- 
gin to form coherent lines, the activation energy increases 
so that V Pt s(l/T) curve downward. However, the fact 
that V p remains greater than V s (below T* , V p /V s rapidly 
increases with decreasing T) indicates that the vortices 
do not extend through the whole sample. 

An important fact for understanding the nature of dis- 
sipation below T* is that the activation energy decreases 
with increasing current and eventually becomes smaller 
than that for 2D pancakes. Figure [| displays the current 
dependence of U^ D = E p at T < T* for H = 0.2 T. The 
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FIG. 3. Current dependence of the activation energy Up D 
determined as dlnV p /dhi(l/T) for T < T* in a field 
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0.2 T. In the limit of small current, Up D « (7, 
shown by the arrow. The value of [/ ps 45 if is also in- 
dicated by the double arrow. Inset: Field dependence of 
the activation energy U% D determined as din V s /dln(l/T) for 
T < T* . The solid line corresponds to H~ ' 2 dependence. 

double arrow indicates the value of the current indepen- 
dent U 2D = E p>s at T > T* and for the same magnetic 
field. The fact that U% D < U 2D for I > 0.2 mA indicates 
that the vortices remain coupled even at large currents 
which allows a new channel of relaxation to emerge and 
become dominant due to smaller activation energy. The 
inset to Fig. is a log-log plot of U^ D vs H, obtained 
from the secondary voltage V s (limit of very small cur- 
rent). The values of U^ D follow an iJ -1 / 2 dependence. 
This field dependence of Ug D is consistent with a 3-D 
plastic creep model, where Uf D is estimated as the en- 
ergy needed for the formation of a double kink over the 
Peierls barrier |14|,[l5| . These results show that the dissi- 
pation at T < T* is determined by two parallel processes: 
thermally activated motion of correlated vortices (domi- 
nant at low currents) with the activation energy greater 
than that for 2D pancake vortices, and plastic motion 
of dislocations (dominant at higher currents) with the 
activation energy smaller than that for a 2D vortex. 

At lower temperatures, the primary voltage V p be- 
comes temperature independent for T < T q (H) and 
scales approximately with the field [Fig. [j](a)] while V s re- 
mains thermally activated, but with a noticeably smaller 
activation energy (Fig. Bj. 

The high normal state resistivity of this specimen may 
be the reason that quantum creep begins to dominate 
the thermally-activated process at relatively high tem- 
peratures T q ~ 5 K. A previous study has shown 
that the normal state resistivity p n (T) of this crystal 
(revealed by the suppression of superconductivity with 



a large magnetic field) is insulating, similar to that of 
PrBa 2 Cu 3 7 ^ s , so that p n (T) -> oo at T -> fil|. The 
reduced dissipation in the normal core of a vortex due to 
a large normal-state resistivity increases the mobility of 
the vortices and facilitates tunneling 113 ]. 

The succession of the coupling transition by quan- 
tum creep points towards the model of flux flow shown 
schematically in the inset to Fig. 0(b) . Below T* the vor- 
tices begin to form coherent lines and the dissipation be- 
comes non-ohmic. The transport current is mostly con- 
fined within a narrow layer below the current contacts. 
The rest of the sample remains relatively undisturbed. 
These two regions of the sample are uncoupled and, inside 
each of them, the vortices are coherent over a distance 
comparable to the thickness of the respective layer. In 
this scenario, the large resistive anisotropy in the mixed 
state results from the loss of the phase coherence only 
between two macroscopic regions of the sample, not be- 
tween all microscopic layers (such as C11O2 bilayers) as 
in the 2D phase at T > T*. The fact that at T < T q we 
observe a T-independent V p and a thermally activated 
V s indicates that the thickness of the current-carrying 
layer is just a few unit cells (which makes it similar to 
ultra thin films and multilayers, the only other systems 
in which quantum creep was observed in transport mea- 
surements). The small length of vortices makes tunneling 
a dominant mechanism even at relatively high tempera- 
tures (T k 5 K); the vortices in the lower layer (which 
is the bulk of the sample) do not tunnel because of the 
much longer correlation length and, therefore, V s remains 
thermally activated. 

This scenario raises the possibility that the trans- 
port current undergoes self-channeling at T < T* . 
In the normal and mixed states above T* , the 
Y\- x Pr x BaiCuj,0'j-t, crystals are not very anisotropic; 
for the crystal studied (length L w 1 mm and thick- 
ness D w 0.015 mm), the ratio V p /V s - 2 at T > T* , 
so that the transport current fills the cross-section fairly 
uniformly. On the other hand, the transition to quantum 
creep, which is extremely sensitive to the length of the 
tunneling segments, indicates that near T q , the current- 
carrying volume collapses into a layer just a few unit cells 
thick. A possible reason may be the non-linear as well 
as nonlocal relationship between the current density and 
electric field in the mixed state below T* . The non-linear 
resistive anisotropy may result in instability of the cur- 
rent distribution with respect to formation of a very thin 
current-carrying channel. 

In order to confirm that the T-independent V p 
[Fig. |l|(a)] is the result of quantum flux creep, we per- 
formed magnetic relaxation measurements on a similar 
crystal. At higher temperatures the relaxation of the 
magnetic moment proceeds as a sequence of uncorrelated 
microscopic steps, each involving thermal activation over 
a certain energy barrier U. The decay time r^ during 
which the induced moment loses a substantial fraction of 



its initial value can be defined as: 

(U(H,T) 
Td = r esc exp <^ — 



(1) 



where the Boltzman factor reflects the degree of avail- 
ability of energy U required for an average elementary 
step and is essentially independent of the physics of the 
relaxation process, while the pre-exponential factor r esc 
is a measure of how rapidly the relaxation would proceed, 
had it not been limited by unavailability of thermal en- 
ergy. We call r esc the escape time, to distinguish it from 
the microscopic attempt time r which characterizes the 
period of vibration of the vortex inside a pinning well. 
The escape time may depend on the sample size, mag- 
netic field, and temperature. Factorization of Td given 
by Eq. (|l|) is meaningful as long as the Boltzman factor 
exp(C//T) 3> 1, so that it dominates the T and H depen- 
dence of T rf . With U(H,T) = U (H)(1 - T/T cr ), where 
T cr is the temperature at which the activation energy 
vanishes, the decay time becomes: 



Td 



.exp 



'•■<"»(?-£ 



(2) 



We emphasize that Eq. (|J) is more general than any par- 
ticular dynamic model of the relaxation process driven 
by fluctuations. 

Within the decade of time 10 3 — 10 4 s, the magnetic 
relaxation curves are well fitted to M{ rr — a — bln(t/to), 
where to is an arbitrary unit of time. Due to the slowness 
of the relaxation, the decay time cannot be determined 
directly by monitoring the induced moment until it loses 
a substantial fraction of its initial value. An alternative 
method, which we will use, is to estimate the decay time 
by extrapolating the Mi rr (t) dependence to Mi rr (rd) = 
0, which gives: 



Td = t exp{a/b}. 



(3) 



With this definition, Td is universal and does not depend 
on the choice of to so that we can compare the experi- 
mental Td with Eq. pi. Figure 111(b) shows Td calculated ac- 
cording to Eq. (H) and plotted vs 1/T for different values 
of H . At higher temperatures the data display an Arrhc- 
nius dependence with a slope dhiTd/d{l/T) decreasing 
with increasing H . The trend is consistent with the field 
dependence of the activation energy in transport mea- 
surements, Fig. 0(a). 

According to Eq. |2J, the pre-exponential factor T esc can 
be determined by the extrapolation of the Arrhenius de- 
pendence of Td(T) to T = T cr . Indeed, the straight 
lines extrapolating the activated dependence converge at 
T cr w 19 K s=s T c and Td — T esc s=s 1 s. The value of 
T esc is consistent with the time it takes for a vortex to 
diffuse from the bulk to the outer edge of the sample: 
T esc ~ R 2 jV v ~ R 2 /u; a £ 2 , where R is the characteristic 



size of the sample in the direction of diffusion, and T> v is 
the diffusion coefficient determined by the attempt fre- 
quency io a and the average hopping distance t a . With 
T esc ~ 1 s and R 2 ~ 10~ 2 — 10 -3 cm 2 (for the crystal 
we used), V v - 1(T 2 - 1CT 3 cm 2 /s. This value of V v 
is consistent with an elementary step of the order of the 
correlation length £ a ~ 100 A and uj a ~ 10 11 - 10 12 s _1 . 

At lower temperatures, the decay time saturates at a 
roughly T- and H- independent level [Fig. 0(b)]. The val- 
ues of T q (H) from transport [(Fig. ffl(a)] and magnetic re- 
laxation [Fig. Mb)] measurements are very close in spite 
of a very large difference in the currents involved in these 
measurements. The fact that the transition to a temper- 
ature independent dissipation takes place in both trans- 
port and magnetic relaxation processes, and at approx- 
imately the same temperature in a given field indicates 
that both phenomena have the same origin. The Euclid- 
ian action Se can be estimated as S^/ft — \n(Td/T esc ) ss 
25, which is comparable, but somewhat smaller than pre- 
viously reported values for other systems j| . 

In summary, we have observed quantum creep in 
underdoped Y\- x Pr x Ba<iCu?,0>i-b crystals using both 
transport and magnetic relaxation measurements. The 
transition to quantum creep is preceded by a coupling 
transition T* which leads to non-ohmic dissipation. The 
transport data indicate that, at T < T* , most of the cur- 
rent is confined to a layer just a few unit cells thick below 
the current contacts. This current-carrying layer is de- 
coupled from the rest of the crystal, where the vortices 
are mostly undisturbed by the current and are coherent 
over a macroscopically long distance. The transition from 
thermally activated to temperature independent dissipa- 
tion takes place only in the current-carrying layer. The 
relatively large probability of tunneling which makes pos- 
sible the observation of this crossover at T m 5 K is due 
to the very short length of the tunneling segments. 
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